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^ ■ ABSTRACT 

o : 

^ . We have mapped, for the first time, the full velocity extent of the water maser 

emission in NGC 3079. The largely north-south distribution of emission, aligned with 
a kpc-scale molecular disk, and the segregation of blue- and red-shifted emission on the 
sky are suggestive of a nearly edge-on molecular disk on pc-scales. Positions and line- 
' of-sight velocities of blue- and red-shifted maser emission are consistent with a central 

Q>^ ' mass of ~ 2 X 10 Mq enclosed within a radius of ~ 0.4 pc. The corresponding mean 

'nI" ■ mass density of 10^'® Mq pc~^ is suggestive of a central black hole, which is consistent 

QQ . with the detection of hard X-ray excess (20 — 100 keV) and an Fe Ka line from the 

, nucleus. Because the rotation curve traced by the maser emission is flat, the mass 

! of the pc-scale disk is significant with respect to the central mass. Since the velocity 

"fH \ dispersion of the maser features does not decrease with radius and constitutes a large 

Qh| fraction of the orbital velocity, the disk is probably thick and flared. The rotation curve 

O ' and the physical conditions necessary to support maser emission imply a Toomre Q- 

' parameter that is ^ 1. Thus, the disk is most likely clumpy, and we argue that it is 

d • probably forming stars. Overall, the accretion disk in NGC 3079 stands in contrast to 

J> . the compact, thin, warped, differentially rotating disk in the archetypal maser galaxy 

; NGC 4258. 

^ I We have also mapped radio continuum emission in the vicinity of the disk and 

■ " " ' identify a new, time- variable, non-thermal component (E) that is not collinear with the 

previously imaged putative jet. Based on the large luminosity and the unusually steep 
spectrum (a < —2.1), we exclude a radio supernova as the progenitor of E. However, 
because its spectrum is consistent with an aging electron energy distribution, E might be 
a rapidly cooling remnant, which may indicate that the jet axis wobbles. Alternatively, 
considering its location, the component might mark a shock in a wide-angle outflow that 
is interacting with a dense ambient medium. In this context, masers at high latitudes 
above the disk, mapped in this and previous studies, may be tracing an inward extension 
of the kpc-scale bipolar wide-angle outflow previously observed along the galactic minor 
axis. 



Subject headings: galaxies: active — galaxies: individual (NGC3079) — galaxies: 
Seyfert — galaxies: starburst — ISM: jets and outflows — masers 



2 



1. Introduction 

NGC 3079 is a relatively nearby (vgys = 1125 ± 6kms~^), edge-on spiral galaxy classified 
as Seyfert 2 (Ford et al. 1986; Ho, Filippenko, & Sargent 1997) or a LINER (Heckman 1980). 
Observations of large-scale structure in this galaxy indicate the presence of a powerful outflow that 
appears to originate within the nucleus. Along the minor axis of the galaxy, at a position angle 
(P. A.) of ~ 80°, a bipolar super-wind inflates a kpc-scale superbubble visible in radio continuum 
emission (Duric & Seaquist 1988; Baan Sz Irwin 1995; Irwin &: Saikia 2003), [NIIJ+Hct emission 
(Veilleux et al. 1994), and soft X-ray continuum emission (Cecil, Bland-Hawthorn, & Veilleux 2002). 
The apex of the superbubble is embedded in a dense molecular kpc-scale disk traced by CO (Koda 
et al. 2002), HCN, and HCO+ emission (Kohno et al. 2001). The circumnuclear molecular disk is 
aligned with, and rotates in the same sense as, the much larger galaxy disk such that the north 
side is approaching and the near side is projected to the west of its major axis. 

There is compelling evidence that NGC 3079 contains an active galactic nucleus (AGN). The 
X-ray spectrum exhibits a continuum excess in the 20 — 100 keV band and a 6.4 keV Fe Ka line, 
which are both believed to be unambiguous indicators of nuclear activity (lyomoto et al. 2001; 
Cecil et al. 2002). Both soft and hard X-ray data are suggestive of unusually high hydrogen column 
density towards the nucleus, Nh ~ 10^^ cm~^ (lyomoto et al. 2001), which, in the context of the 
AGN unified model, is indicative of an almost edge-on obscuring structure along the line of sight 
to the nucleus (Lawrence & Elvis 1982; Antomicci 1993). The existence of nuclear activity is also 
indirectly supported by radio observations of the nucleus. Continuum images of the inner-pc region 
obtained with Very Long Baseline Interferometry (VLBI) reveal a linear distribution of compact 
knots consistent with a jet (Irwin &; Seaquist 1988). Furthermore, the galaxy harbors one of the 
most luminous (i.e., assuming isotropic emission of radiation) water masers known (Henkel et al. 
1984; Haschick &; Baan 1985), and such maser emission generally is associated with AGN activity 
(Braatz, Wilson, & Henkel 1997). In fact, the isotropic luminosity of the brightest, spatially 
unresolved spectral feature in NGC 3079 (~ 130 Lq) is comparable to luminosities of other AGN 
related masers such as NGC 4258 and NGC 1068 but at least two orders of magnitude greater than 
the isotropic luminosity corresponding to the total integrated maser flux from starburst galaxies 
such as M82 (Baudry, Brouillet, & Henkel 1994) and NGC 253 (Ho et al. 1987) and from the W49N 
maser, the most luminous maser in our galaxy associated with a region of intense star formation 
(Gwinn 1994). 

The nucleus of NGC 3079 may also host a starburst; the most convincing evidence comes 
from infrared (IR) and far-infrared (FIR) observations. Lawrence et al. (1985) found color excess 
and an extended lO^m nuclear source indicative of large amounts of dust in the nuclear region. 
The IRAS 60 to 100 /im flux ratio is ~ 0.5, which is suggestive of warm dust on kpc scales and 
qualifies NGC 3079 as a starburst galaxy (Soifer et al. 1989; Braine et al. 1997). FIR observations 
of the galaxy with ISOPHOT reveal a powerful ~ 3.0 x 10^° Lq point source of FWHM < 4.5 kpc 
coincident with the nucleus (Klaas &; Walker 2002). The SED of this point source is indicative 
of 20 - 32 K dust temperature. The luminosity of 3.5 x 10^° from IRAS (Soifer et al. 1987) 
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or ~ 3.0 X 10^° Lq from ISOPHOT (Klaas k Walker 2002) yields a star formation rate of ~ 10 

Mq yr~-^ on kpc scales (Veilleux et al. 1994), under the assumption that star formation alone 
contributes to the infrared flux. The presence of nuclear star formation may also be inferred from 
large quantities of molecular gas associated with the nucleus. The observations of CO emission 
indicate a dynamical mass of 7 x l(f Mq enclosed within a radius of 76 pc (Koda et al. 2002), 
as well as an unusually large molecular surface gas density relative to other normal, starburst, or 
active galaxies of ~ 7200 Mq pc~^ (Planesas, Colina, Sz Perez-Olea 1997; Sakamoto et al. 1999; 
Koda et al. 2002). Although sufficient to fuel a starburst, such a high concentration of gas does 
not necessarily imply intense nuclear star formation (Jogce 2001). However, the dynamical and gas 
masses deduced from CO observations by Koda et al. (2002) are suggestive of a Toomre-Q parameter 
that is less than unity within 150 pc of the nucleus and for isothermal sound speeds < 50kms~^, 
consistent with the presence of molecular gas. Thus, the unusually high concentration of gas in 
the nucleus of NGC 3079 is subject to gravitational instabilities and might be actively forming 
stars. Detailed studies of the kpc-scale superbubble arc also indicative of the nuclear starburst. 
The inferred starburst age (> lO'^ years) and gas depletion time (> 10^ years) are longer than the 
dynamical age of the superbubble {^^ 10^ years), which is consistent with a starburst driving the 
bubble expansion (Veilleux et al. 1994). Furthermore, based on (1) the similarity of the superwind 
in NGC 3079 to that in prominent starburst galaxies like M82 and NGC 253 and (2) the general 
lack of correlation between the properties of outflow and nuclear activity in a sample of galaxies, 
Strickland et al. (2004a) argue that supernova feedback is responsible for the superbubble in NGC 
3079. The correlation between diffuse halo X-ray luminosity per unit stellar mass and the infrared 
luminosity per unit mass in a sample of seven starburst and three "normal" galaxies (see Fig.4 in 
Strickland et al. 2004a) is supportive of the starburst hypothesis, although that correlation does 
not imply causation particularly in the context of an alleged starburst-AGN connection. 

Arguments have been made against the existence of a nuclear starburst in NGC 3079. The 
ratio of FIR luminosity to the molecular gas mass (i.e., the star formation efficiency) computed 
on kpc scales is Lfir/Mh2 ^ 13, which is greater than the value of ~ 4 for the Milky Way and 
1 < Lpjr/Mh2 < 7 for Giant Molecular Clouds, but below the average value of ~ 25 observed for 
starburst galaxies (Young et al. 1986; Planesas et al. 1997). In fact, small-beam (6") 10 //m flux 
of NGC 3079 measured by Lawrence et al. (1985) yields Lfir/Mh2 of only 4.7 on smaller scales, 
which is suggestive of low star formation efficiency in the nuclear region (Hawarden et al. 1995). 
Furthermore, since the 10 jim. flux in a 5.5" aperture is only 4% of the larger beam color corrected 
IRAS flux, the infrared emission is not centrally concentrated, contrary to expectation for a nuclear 
starburst (Devereux 1987). The nuclear region of NGC 3079 (< 500 pc) is characterized by a flat 
spectrum between 10 and 20 /xm in contrast to what is observed for starburst galaxies (Lawrence 
et al. 1985; Hawarden et al. 1986; Klaas &; Walker 2002). According to Hawarden et al. (1995), low 
extinction towards the parts of the nucleus not obscured by the molecular disk, the lack of a 10 /xm 
compact source coincident with the relatively well-defined apex of the superbubble, as well as the 
unusually high H2/Br7 ratio all argue against the starburst hypothesis. However, the validity of 
these arguments may be in doubt because of extremely heavy obscuration in the nuclear region due 
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to dust lanes and the massive kpc-scale molecular disk. Indeed, asymmetry in the K-band image 
of the nucleus implies considerable extinction even at 2.2 /xm (Israel et al. 1998). 

If the nuclear starburst hypothesis is accepted, then it is not unreasonable to suppose the 
coexistence of star formation and nuclear activity in the inner parscc of NGC 3079, although 
the particular relationship between the two remains uncertain. It has been suggested that star 
formation and nuclear activity are likely to coexist in galactic nuclei, either because they are 
coupled through evolutionary mechanisms, or simply because they both depend on gas inflow and 
accretion (Cid Fernandes et al. 2001, and references therein). The existence of star formation and 
of young stellar populations with ages of (7—20) x 10^ years in the circumnuclear tori of radii < 1 pc 
has been proposed to explain the UV excess of many Seyfert 2 galaxies (Tcrlevich 1996). Our galaxy 
contains a population of young, hot, bright, and massive stars within the inner-pc of the central 
black hole (Genzel et al. 1996, 2003), which might have formed in-situ in a preexisting accretion 
disk (Milosavljevic &; Loeb 2004). A young stellar population with an age of ~ 7x 10^ years has been 
found in the central 12 pc of the Circinus active galaxy (Maiolino et al. 1998), and young massive 
stars with ages < 10^ years have been unambiguously detected within ~ 200 pc of 30% to 50% of 
Seyfert 2 nuclei surveyed (Mas-Hcssc ct al. 1994, 1995; Cid Fernandes et al. 2001, and references 
therein). Because of the likely coexistence of star formation and nuclear activity in NGC 3079, the 
proximity of the galaxy, the presence of maser emission, and apparent jet activity, NGC 3079 is 
an important case to study with regard to understanding connections between star formation and 
nuclear activity. 



1.1. Previous Imaging Studies of the Nucleus 
1.1.1. Continuum 

VLBI continuum images of NGC 3079 (refer to Fig.l for a schematic) reveal several coUincar 
components. A, B, C, and D, extending over ~ 2.1 pc (~ 25 mas) at P.A.~ 125° (Irwin & Seaquist 
1988; Trotter et al. 1998; Sawada-Satoh et al. 2000, 2002). A line connecting the continuum 
components is significantly misaligned with respect to the superbubble axis but approximately 
aligned with the southern edge of the biibblc (Cecil et al. 2001). Components A and B are separating 
at about 0.13 c, which strongly indicates the presence of a jet (Trotter et al. 1998; Sawada-Satoh 
et al. 2000, 2002), while the dominant component B is stationary with respect to the molecular gas 
supporting the maser emission to within a transverse speed of 0.02c (Sawada-Satoh et al. 2000). 
These observations led Sawada-Satoh et al. (2000) to identify component B with the central engine. 
However, this interpretation is not unique since, even if component B were a radio jet core, an offset 
of the core from the central engine could be significant (e.g., NGC 4258; Herrnstein et al. 1997). 
Component B would also remain stationary with respect to the maser if B were shock excited 
emission, as would occur where a hot ionized supersonic flow encounters a dense ambient medium 
(e.g., Phillips & Mutel 1982; O'Dea, Baum, & Stanghellini 1991; Aurass, Vrsnak, & Mann 2002, 
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Kellermann et al. 2004). In support of the latter hypothesis, Trotter et al. (1998) found that 
components A and B have similar spectra, exhibit spectral turnovers between 5 and 22 GHz, and 
consist of multiple components or extensions. Since its spectrum is not flat, B is most likely not a 
jet core. Instead, Trotter et al. (1998) suggested that both components are lower redshift and lower 
luminosity counterparts of compact symmetric objects whose flux densities have been observed to 
peak at gigahertz frequencies and which are ascribed to the interaction between jets and associated 
high-density nuclear gas. 

The putative jet coexists with additional non-thermal sources of continuum emission, whose 
nature remains ambiguous. Component E, identified by Kondratko, Greenhill, k, Moran (2000) , and 
component F, identified by Middelberg et al. (2003), both lie ~ 1.5 and ~ 3.7pc approximately east 
of component B, respectively. Middelberg et al. (2003) found that the separation of components E 
and B at 5 GHz systematically decreases at a rate of 0.07±0.01 c (i.e., roughly one beam width over 
~ 3 years). These authors suggested that E is stationary and that component B is moving toward 
E. However, this proposition seems unlikely since B appears to be stationary with respect to the 
extended molecular gas structure underlying the maser emission. A more realistic scenario is that 
component E is a pattern in a flow that only appears to move in the direction of B as its physical 
characteristics change. Such apparently inward motions have been observed in several quasar and 
AGN jets (Kellermann et al. 2004). 



1.1.2. Maser 

In previous observations (Fig. 1), the water maser in NGC 3079 was found to be distributed in 
a disordered linear structure along P.A.~ —10°, aligned with the larger circumnuclear CO disk at 
P.A.~ -11° (Trotter et al. 1998; Koda et al. 2002). Unlike the archetypal water megamaser NGC 
4258, the maser in NGC 3079 displays significant structure orthogonal to the general elongation 
and large velocity dispersion (~ 30knis^-'^) in relatively compact areas on the sky (~ 0.1 pc). 
Although no simple rotation law fits the data well, the velocity distribution of the maser emission 
is consistent with a binding mass of ~ 10^ Mq and indicates that the pc-scale accretion disk rotates 
in the same sense as the kpc-scale molecular disk. The maser predominantly traces the approaching 
side of the disk and thus reveals considerable asymmetry with the blue-shifted emission significantly 
dominating the detected flux density. In particular, the maser emission extends roughly from 890 
to 1190 km s~^ with most emission lying between 930 and 1060 km s~^. In addition to the main 
distribution of the maser emission at P.A.^' —10°, Trotter et al. (1998) found two maser features 
that define a second axis at P.A.~ 128°. Since one of these two maser features coincides on the 
sky with component C and the other is in close proximity to component B, they may constitute a 
second population of water masers in which amplification of the continuum components might play 
a role. 

The difficulty in interpreting the maser data has resulted in two very different models for the 
nuclear region in NGC 3079 (Fig. 1). Relying on the apparent alignment of the distribution of 
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the maser emission with the kpc-scale molecular structure, Trotter et al. (1998) proposed a largely 
north-south, relatively thin, turbulent, and highly-inclined accretion disk for the geometry of the 
nuclear region (Fig. 1). In this model, the dynamical center is located between components A and 
B, while the maser emission occurs on the surface of the disk in shocks stimulated by the interaction 
of a nuclear wind with the disk. The asymmetry in blue- and red-shifted flux density is attributed 
to free-free absorption in the ionized wind, related to the off-axis outflow as traced by components 
A, B, and C. 

The model proposed by Trotter et al. (1998) differs significantly from that proposed by Sawada- 
Satoh et al. (2000). According to these authors, the observed velocity drift of the maser feature 
at 1190 km s~^ is much larger than the velocity drift of the other maser features and is due to 
centripetal acceleration. In a rotating system, only the maser features nearest in projection to 
the rotation axis are expected to show large velocity drifts. Consequently, the rotation axis would 
have to pass close to the 1190 km feature and through component B, the proposed location 
of the dynamical center. In this model (Fig. 1), the masers are embedded in a torus of parsec 
thickness misaligned with respect to the kpc-scale molecular disk by ~ 110°. It is not clear why the 
embedded maser emission should be distributed along the observed locus. However, the outflow 
axis is misaligned with respect to the minor axis of the pc-scalc rotating structure significantly 
more in the Trotter et al. (1998) model than in the Sawada-Satoh et al. (2000) model. 



1.2. Motivation 

The Trotter et al. (1998) and Sawada-Satoh et al. (2000) models for the nuclear region of 
NGC 3079 differ primarily in the location of the dynamical center as well as in the thickness and 
orientation of the rotating structure. It is thus not clear what is the correct geometric model for the 
inner parsec, what stimulates the maser emission, and why it predominantly traces the regions of 
the molecular disk approaching the observer. Ambiguity also remains about the connection of the 
continuum components and the disk and the central engine, about the relationship, if any, between 
the nuclear activity and the starburst, and about the connection between the pc and kpc-scale 
outflows in this galaxy. Our study of NGC 3079 is motivated by these outstanding questions and 
by the disagreements between the two existing models. 

In an effort to better understand the nuclear region in NGC 3079, we attempted to fit a geomet- 
ric and dynamical model of an inclined, thick, turbulent disk to the maser positions and velocities 
reported by Trotter et al. (1998). The model not only reproduced the known emission at velocities 
930 — 1060 km s~^ but also predicted the existence of emission at velocities 1190 — 1320 km s~^ 
located roughly south from the known blue-shifted emission. No previous VLBI observation had 
explored the velocities above 1222 km s^^. However, Nakai et al. (1995) marginally detected a 
0.08 Jy peak at 1267 km s~^ using the 45-m telescope of the Nobeyama Radio Observatory (Fig. 2). 
Motivated by this detection, we imaged with VLBI the water maser emission in NGC 3079 in order 
to determine whether the red-shifted emission lay south of the blue-shifted emission as predicted 
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Fig. 1. — Previously proposed models for the nuclear region of NGC 3079. In the model in the right 
panel, reproduced from Trotter et al. (1998), the maser emission occurs in a thin, almost edge-on, 
and turbulent disk aligned with, and rotating in the same sense as, the kpc-scale molecular disk. 
In the model shown in the left panel, adapted from Sawada-Satoh et al. (2000), the maser emission 
takes place in a very thick, edge-on torus misaligned significantly with respect to the kpc-scale 
molecular disk. In both models, the outflow takes the form of a collimated jet not well aligned with 
the axis of the rotating structure. 



by our model. Recently, more red-shifted features have been detected with the Effelsberg 100-m 
antenna (Hagiwara et al. 2002). In this study, we present the first VLBI map of the red-shifted 
side of the accretion disk in NGC 3079, and thereby hope to clarify the geometry of the inner-pc. 

Data calibration and reduction techniques are discussed in Section 2. Spectral-line and contin- 
uum images of the inner-pc region are presented in Section 3. In Section 4 we interpret the observed 
kinematics of the maser in terms of a thick, flared, disordered, self-gravitating, and clumpy pc-scale 
disk, which orbits a central mass of ~ 2 x 10^ Mq enclosed within 0.4 pc. In this work we assume 
that NGC 3079 is at a distance of 17.3 ± 1.5 Mpc (Tully, Shaya, & Pierce 1992), so that Imas 
corresponds to 0.084 pc. All velocities have been computed in accordance with the radio definition 
of Doppler shift and in the local standard of rest (LSR) frame. 
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Fig. 2. — First possible detection of red-shifted water maser emission in NGC 3079. Spectra 
reproduced from Nakai et al. (1995). 



2. Observations and Calibration 

NGC 3079 was observed in spectral-line mode with the Very Long Baseline Array (VLBA), 
augmented by the phased VLA and the Effelsberg (EE) 100-m antenna, for approximately 12 hours 
on 2001 March 23. The source was observed with 4 x 16 MHz intermediate-frequency (IF) bands, 
which covered a velocity range of 726 to 1491 km s^^. The 5, 8, and 15 GHz radio continuum 
observations were made on 1996 December 2 with the VLBA and the VLA phased array. The source 
was observed with a 32 MHz bandwidth in each polarization and for approximately three hours at 
each frequency. Both spectral and continuum data were processed with the VLBA correlator at 
the NRAO.^ 

The data were reduced and calibrated with standard continuum and spectral-line VLBI tech- 
niques in AIPS^. The amplitude calibration for the two experiments included corrections for at- 
mospheric opacity. Antenna gain curves and measurements of system temperature were used to 
calibrate amplitude data for each of the VLBA stations. Amplitude calibration for the VLA was 
based on VLA scans of 3C 286, for which we adopted a 22 GHz flux density of 2.55 Jy, and on boot- 
strapping of flux densities for VLBI calibrators. To minimize systematic errors due to atmospheric 



^The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooperative 
agreement with the National Science Foundation 

^Astronomical Image Processing System 
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opacity effects, only the VLA scans of VLBI calibrators obtained near the elevation of 3C 286 were 

used in bootstrapping. The data in each polarization were corrected for the effect of time variable 
source parallactic angle, and the station positions were corrected for motions due to plate tectonics 
by utilizing the US Naval Reference Frame solutions 9810 and 9513. 

In the 5, 8, and 15 GHz data, the residual delays and fringe rates due to the troposphere 
and clock uncertainties were removed via observations of calibrators 4C 39.25, 3C273, OJ287, 
and OQ208. The final images of the nuclear region of NGC 3079 were obtained through ap- 
plication of self-calibration. In the spectral data, the complex bandpass shape and any elec- 
tronic phase difference between the IF bands were removed by observation of calibrators 4C 39.25, 
1150-1-812, 1308-1-326, and 3C345. The position of the maser was corrected by a fringe rate anal- 
ysis of 10 spectral channels including the brightest spectral feature at 956kms~^ with the result 
aaooo = 10''01™57!802± 0.001, ^2000 = 55°40'47'.'26±0.01. Phase and amplitude fluctuations due to 
troposphere and clock uncertainties were removed by self-calibration of the brightest maser feature 
at 956kms~^ and application of the resulting solution to all spectral channels in all IF bands. After 
calibrating and imaging the spectral-line data set, the relative positions of the maser features were 
obtained by fitting two-dimensional elliptical Gaussians to the distribution of the maser emission 
in each spectral channel. We obtained a continuum image at 22 GHz by averaging all line-free 
frequency channels. Table 1 lists beam dimensions and noise levels at each frequency. Table 2 lists 
the properties of all spectral-line observations of the water maser with the VLBA. The experiment 
presented in this study images the widest velocity range of emission among all VLBA observations 
made to-date. 

Table 1: Half-power beam dimensions and noise levels. 



Frequency 


Major 


Minor 


P.A. 


RMS 


(GHz) 


(mas) 


(mas) 


n 


(mJy beam~^) 


5 


3.6 


2.4 


14 


0.060 


8 


2.0 


1.3 


29 


0.055 


15 


2.0 


1.3 


-89 


0.11 


22 


0.30 


0.26 


-43 


0.14*^ , - 2.3'' 



"continuum 

''spectral-line for a channel width of 0.42 km s~^. 

To estimate the fraction of the power imaged by the interferometer, we acquired a single- 
dish spectrum of NGC 3079 (see Fig.3) by position-switching in a single polarization mode with 
the NASA Deep Space Network 70-m antenna in Robledo, Spain, on 2003 February 4, using an 
observing setup similar to that described in Greenhill et al. (2003) . After integrating for a total of 
2.5 hours, we achieved an rms noise level in the spectrum of 11 mJy for 1.3kms~^ channel spacing. 
The antenna gain curve was obtained in a single track by measurement of the elevation dependence 
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Table 2: VLBA sLudies of llie water masor in NGC 3079. 





Trotter et al. 1998 


Sawada-Satoh et al. 2000 


This study 


Date 


1995 Jan 9 


1996 Oct 20 


2001 Mar 23 


V Range*^ 


880 - 1204 


590 - 1008 


726 - 1491 


AV^ 


0.21 


0.42 


0.42 




0.22 


0.54 


0.14 


ad 


2.8 


6.6 


2.3 



"Total velocity coverage in kms~ 
'Channel spacing in kms~^ 

'^The rms noise level in the pseudo continuum image in Janskys 
''The rms noise level in the spectral-line image in Janskys 

of 1308+326 antenna temperature, corrected for atmospheric opacity estimated from a tipping 
scan. The antenna efficiency was estimated to be 0.43 ± 0.01 based on several measurements of 

antenna temperature for 3C286. From the rms deviation of our gain measurements about the best 
fit polynomial and the formal uncertainty in the efficiency, we estimate that the gain calibration of 
the antenna is accurate to within 10%. 

3. Results 
3.1. Maser 

In agreement with previous observations of NGC 3079, we found that the maser exhibits a 
conspicuous asymmetry whereby the blue-shifted emission significantly dominates the detected flux 
density (Fig. 3). The spectrum of the imaged power agrees satisfactorily with the spectrum detected 
with the Effelsberg 100-m antenna 10 days prior to the VLBI observation (see Fig. 2c in Hagiwara 
et al. 2002). The observation with the NASA Deep Space Network 70-m antenna in Robledo, Spain, 
conducted 1.9 years later showed significant flux variability throughout the spectrum except for the 
strongest feature (Fig. 3). Such variations, however, are not unusual since single-dish monitoring 
of the maser has revealed substantial flux variability on timescales of years (Nakai et al. 1995; Baan 
k Haschick 1996; Hagiwara et al. 2002). 

As with the previous investigations (Trotter et al. 1998; Sawada-Satoh et al. 2000), we found 
that the maser emission is distributed in a disordered linear structure at P.A.~ —10° (Fig. 4). 
Furthermore, the distribution of emission on the sky is clearly segregated by velocity with emission 
blue- and red-shifted with respect to the systemic velocity located in the northern and southern 
section of the image, respectively. To determine temporal changes in the maser distribution, we 
compared the maser positions obtained here with those of Trotter et al. (1998). Both the 1995 
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Fig. 3. — Imaged-power and total-power spectra of the maser in NGC 3079. The total-power 
spectrum (gray) was obtained by position-switching with a 70-m Deep Space Network antenna in 
Robledo, Spain on 2003 February 4. The spectrum of the imaged power (black line) is based on the 
VLBI observation conducted on 2001 March 23 and was computed by summing the flux densities 
of the individual maser spots. Left inset: Imaged-power spectrum (black line) of the brightest 
maser feature at 956kms~^ and a Gaussian fit to the main peak (gray) corresponding to an 
isotropic luminosity of 131 Lq. Right inset: Total-power and imaged-power spectra of the red- 
shifted emission with the vertical scale expanded and the total-power spectrum shifted upward to 
facilitate comparison. The vertical arrow indicates the systemic velocity of the galaxy (1125 km s~^). 



January 9 (Trotter et al. 1998) and the 2001 March 23 (this study) epochs were calibrated by 
referencing the phases of all spectral channels to the brightest maser feature at ~ 956kms~^ and 
thereby arc most likely registered to the same volume of masing gas. Indeed, three-dimensional 
flux density-weighted correlation of the distribution of the maser emission presented in this study 
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Fig. 4. — Distribution of maser emission (colored circles) and of the 22 GHz continuum (gray 
contours) in the nuclear region of NGC 3079. The color of the maser spots indicates line-of-sight 
velocity in accordance with the bar on the right. The continuum contour levels are —3a, 3a, 5a x 
2^^/-^, where a = 0.14 mJy beam~^ and k = 0, 1, 2, 3, . . .. Contours above 5a are shaded. The gray 
curve encloses the newly imaged emission, while the ellipse in the bottom right corner of the figure 
illustrates the resolution beam. 



with that presented in Trotter et al. (1998) indicates that the two epochs are aligned to within 
a = 3.2 kms^^ in velocity and are registered by the self-calibration technique to within a = 0.5 mas 
on the sky (0.5 mas corresponds to a transverse speed of 0.02 c over 6.25 years and is not an intrinsic 
measurement limit on the maser proper motions). In fact, spectra of nearly half of the maser 
emission clumps imaged in both Trotter et al. (1998) and this study peak within 3.2 kms~^ of each 
other, which is comparable to the a = 2.4kms~^ half width of the brightest maser feature (Fig. 3). 
The comparison of the two epochs (Fig. 5) indicates that (1) the newly imaged red-shifted emission 
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Fig. 5. — Angular distribution of the water maser emission in NGC 3079 obtained on 1995 January 
9 (Trotter et al. 1998) and on 2001 March 23 (this study). The top and bottom panels show 
emission blue- and red-shifted with respect to the systemic velocity, 1125 km s~^, respectively. The 
velocity and angular scales are expanded to facilitate comparison between the two epochs. 



is confined to a region south to south-west from the previously known blue-shifted features, (2) 
several new blue-shifted features have appeared, (3) the maser feature located at (5.9, —10.5) mas 
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appears to have drifted in velocity from 1123 to 1139 km s~ , a change of ~ 16kms~ in 6.25 years, 
and (4) the maser feature located at velocity ~ 940kms~^ and at position (—1.4,0.7) mas appears 
to have moved westward with velocity of ~ 0.02 c (~ 6000 kms^^). The latter two cases might each 
be due to variability of separate and distinct gas clumps rather than physical motions. Note that the 
maser feature that has drifted in velocity between Trotter et al. (1998) and this study is significantly 
displaced on the sky from the ~ 1190 km s~^ maser feature [located at (3.1, —14.9) mas], on whose 
velocity drift the model of Sawada-Satoh et al. (2000) depends. Based on the Trotter ct al. (1998) 
and Hagiwara et al. (2002) studies, we place an upper limit of 0.5kms~^ yr"^ on the velocity drift 
of the ~ 1190 km s~^ feature. 



3.2. Continuum 

Components A and B were detected at all frequencies while component C was detected at 
5, 8, and 15 GHz, but not at 22 GHz (Fig. 6, 7, 8, and 4). In agreement with previous reports, 
components A, B, and C are collinear along the position angle of ~ 126°. To estimate the flux 
density, position, and position uncertainty of each continuum component, we fitted multiple two- 
dimensional Gaussian model components to the emission peaks at each frequency (Table 3). In 
agreement with previous reports (e.g., Sawada-Satoh et al. 2000) and supplementing our data with 
published positions, we found that the systematic motion of component A at 22 GHz is consistent 
with a transverse speed of 0.12 ± 0.02 c along P. A. of ~ 126° ± 5°, which is approximately away 
from component B (Fig.9). Note that the change (typical magnitude of ~ 0.35 mas) from epoch 
to epoch in the position of component B with respect to the reference maser feature does not 
appear to be systematic and is inconsistent with formal measurement uncertainties (Fig. 9). This 
observed scatter might be due to the change in morphology or position of either component B or 
the reference maser feature. However, the time average position of component B is approximately 
stationary with respect to the maser emission. We emphasize that the jitter of component B, like 
the motion of E in the direction of B reported by Middelberg et al. (2003), might not be a real 
but an apparent motion due to changing source characteristics in a flow or across shocks (e.g., 
Keller mann et al. 2004). 

We confirmed the spectral turnover reported by Trotter et al. (1998) of components A and B 
and constrain the frequency of this turnover to the range between 5 and 15 GHz (Fig. 10). To place 
more stringent constraints on the turnover frequencies, we considered two simple spectral models for 
each component. In the first model, we assumed that at lower frequencies the spectra are dominated 
by synchrotron self-absorption (Si, oc z/^/^) intrinsic to the source, while at higher frequencies 
the spectra are determined by optically-thin synchrotron emission (3,^ oc v"') also internal to the 
source. Under these assumptions, the spectra of components A and B can be characterized by 
Sv = av^f^^i — e"*'^" ^''^), where a, 6, and a are parameters of the model. The fiux densities 
measured at three frequencies for each component uniquely determine the three parameters of the 
model, which yielded a = —1.5 and —1.1, and turnover frequencies (i.e., frequency of maximum 
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Fig. 6. — Contours of the 5 GHz continuum drawn at levels —2a, 4cj, 6cj x 2^^'^, where a = 
0.060 mJybeam"^ and k = 0,1,2,3,.... Contours above 6o" are shaded. Compact component 
D detected by Trotter et al. (1998) at 5 GHz with flux density of 4.3 ± 0.5 mJy would be located at 
(38.3, —31.9) mas in the above plot but was not detected at a level of 6ct = 0.36 mJy. The maser 
features shown here as small filled circles were registered to the 5 GHz continuum by aligning the 
positions of component B at 5 and 22 GHz. The uncertainties in this registration are 0.02 and 
0.03 mas in right ascension and declination, respectively. The conversion from the color of the 
maser spots to their line-of-sight velocity is given by the bar on the right. Beam dimensions are 
shown in the bottom right corner. 



flux density) of 9.0 and 9.6 GHz for components A and B, respectively. For the data reported by 
Trotter et al. (1998), we obtained a = —1.5 and —1.3, and turnover frequencies of 9.1 and 9.0 GHz 
for components A and B, respectively. However, if free-free absorption [tu oc z/~^'^) extrinsic to 
the source dominates at lower frequencies instead of synchrotron self-absorption, then the spectra 
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Fig. 7. — Contours of the 8 GHz continuum drawn at levels —3a", 3a, 6a x 2^^/^, where a = 
0.055 mJybeam"^ and k = 0,1,2,3,.... Contours above 6a are shaded. The maser features are 
shown using small filled circles color-coded by velocity in accordance with the bar on the right. The 
registration of the maser to the 8 GHz continuum was achieved by aligning the positions of com- 
ponent B at 8 and 22 GHz. The uncertainty in this registration is 0.01 mas in both right ascension 
and declination. The ellipse in the bottom right corner of the figure illustrates beam dimensions. 



of components A and B can be modelled by S,^ = av'^e"^^ ^ (e.g., Bicknell, Dopita, & O'Dea 
1997; Tingay &: Murphy 2001). For this model, we obtained a = —1.8 and —1.3 as well as turnover 
frequencies of 8.7 and 9.5 GHz for components A and B, respectively. Similarly, the flux densities 
reported by Trotter et al. (1998) yielded a = —1.8 and —1.6 as well as peak locations of 8.8 and 
8.7 GHz. Thus, for both models and for both data sets, spectral turnover frequencies lie in the 
range 8-10 GHz. 

The spectra of components A and B are reminiscent of spectra of gigahertz peaked spectrum 
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Fig. 8. — Contours of the 15 GHz continuum drawn at levels —2a, 4cr, 6a x 2*^/^, where a = 
0.11 mJybeam"^ and k = 0, 1, 2, 3, . . .. Contours above 6a are shaded. Superposed on top of the 
continuum are the maser features registered to the 15 GHz continuum by aligning the positions 
of component B at 15 and 22 GHz. The uncertainty in this registration is 0.01 mas in both right 
ascension and declination directions. The bar on the right indicates the relationship between the 
color and the velocity of the maser features. Beam dimensions are shown in the bottom right corner 
of the figure. 



(GPS) sources, as already noted by Trotter et al. (1998). These sources are thought to arise from 
regions where jets interact with the dense ambient gas (e.g., Bicknell et al. 1997). Their spectra 
peak in the frequency range of 0.1 — 10 GHz with high-frequency spectral indices of —1.3 ^ a < 
—0.5 and low-frequency indices of a ~ 1. Because of the similarity between the two components 
and the GPS sources, we suggest that components A and B arise in regions where a relativistic 
outflow is interacting with a dense ambient medium. Since components A and B consist of multiple 
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Table 3: Characteristics of continuum components at 5, 8, 15, and 22 GHz obtained by two- 
dimensional multiple component Gaussian fits to emission features. 



Frequency 


Component 


5. 


R 




Major 


Minor 


P.A. 


(GHz) 




(mJy) 


(mas) 


D 


(mas) 


(mas) 


(°) 


5 


B 


13.4 ±0.4 


— 


— 


3.63 ± 0.07 


1.62 ± 0.03 


19 ± 1 




A 


4.7 ± 0.3 


25.25 ± 0.05 


126.5 ± 0.1 


3.0 ± 0.1 


1.60 ± 0.07 


20 ± 3 




C 


3.6 ±0.5 


17.8 ±0.1 


123.5 ±0.4 


3.4 ±0.4 


2.7 ±0.3 


109 ± 19 




E 


2.8 ±0.6 


25.9 ±0.2 


85.5 ±0.5 


4.0 ±0.7 


2.7 ±0.4 


32 ±17 


8 


B 


32.4 ±0.2 






2.89 ±0.01 


1.87 ±0.01 


24.2 ±0.3 




A 


11.5 ±0.2 


25.07 ± 0.02 


126.94 ± 0.03 


3.11 ±0.04 


2.31 ± 0.03 


140 ±2 




C 


9.7 ±0.3 


16.81 ±0.03 


124.4 ±0.1 


4.5 ±0.1 


2.58 ± 0.06 


135 ±2 




E 


2.5 ± 0.2 


25.90 ± 0.06 


87.4 ± 0.2 


2.9 ±0.2 


2.5 ±0.2 


175 ± 20 


15 


B 


27.3 ±0.3 






2.71 ±0.02 


1.95 ±0.02 


97.6 ± 0.6 




A 


7.8 ± 0.5 


24.84 ± 0.09 


127.2 ±0.2 


4.6 ±0.3 


2.3 ±0.1 


128 ±3 




C 


3.7 ±0.4 


16.91 ±0.09 


124.6 ± 0.3 


3.0 ± 0.2 


2.0 ±0.2 


130 ±7 


22 


B 


40 ±1 






1.09 ± 0.03 


0.88 ± 0.03 


31 ±6 




A 


28 ±2 


27.14 ± 0.03 


127.96 ±0.06 


1.24 ± 0.06 


1.08 ± 0.05 


24 ±14 



Note. — Angular separation R and position angle <j) are measured with respect to component B at each frequency. 
The 5, 8, and 15 GHz data were acquired on 1996 February 12 while the 22 GHz data were obtained on 2001 March 
23. Component D (S^ = 4.3 ± 0.5 mJy) at position r = 54.5 ± 0.4 mas and (f) = 124.7° ± 0.4° and detected by Trotter 
et al. (1998) was not detected at a level of 6(t = 0.36 mJy in this study. One mas corresponds to 0.084 pc. 

components and do not exhibit a flat spectrum, neither A nor B is likely to be a jet core (cf. 
Sawada-Satoh et al. 2000). 

Trotter et al. (1998) reported a detection at 5 GHz with a signal-to-noise ratio (S/N) of 39 of a 
compact 4.3 lb 0.5 mJy component (D) located ~ 4 pc southeast of B and collinear with components 
A, B, and C. However, we found no evidence above the 0.36 mJy (6(7) level of this component at 
5 GHz. Since the 5 GHz beam area of Trotter ct al. (1998) was 2.4x larger than the 5 GHz beam area 
in our study, the non-detection of D could be due to either resolution efl'ects or source variability. 
Component D was also not detected at 5 GHz by Middelberg et al. (2003). 

The continuum images at 5 and 8 GHz revealed a new component (E) that is not collinear 
with previously known components A, B, C, and D (Fig. 6 and Fig. 7). We detected component 
E with an S/N of 31 and 16 at 5 and 8 GHz, respectively (Kondratko et al. 2000). Middelberg 
et al. (2003) confirmed the detection of E and also reported identification at 1.7 and 2.3 GHz of a 
new component F located roughly 3.7 pc from component B along P.A. of ~ 100°. Since the beam 
size in the Trotter et al. (1998) study was larger than the beam size in our study, component E 
should have been detected by Trotter et al. (1998) with an S/N of at least 17 at 5 GHz if its flux 
density were constant. The failure of Trotter et al. (1998) to detect E indicates that the source 
is most likely time variable. Indeed, comparison of our results with those of Trotter et al. (1998) 
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Fig. 9. — Left: Position of component B with respect to the reference maser feature. The small 
scatter (typical magnitude of ~ 0.35 mas) in the relative position of component B is not consistent 
with formal measurement uncertainties and might be due to morphological or positional changes 
in the reference maser feature or in component B itself. Top right: Position of A with respect to 
the reference maser feature. A fit to the four measured positions indicates that A is moving along 
P. A. ~ 126° lb 5°, effectively away from component B. Bottom right: Temporal motion of A along 
the linear trajectory in the top right panel, consistent with a transverse speed of 0.12 it 0.02 c. The 
1995 January 9 epoch was obtained from Trotter et al. (1998), epochs 1996 October 20, 1997 March 
7, and 1997 October 2 are from Sawada-Satoh et al. (2002), while the 2001 March 23 epoch is from 
this study. For points where error bars are not visible, the errors are smaller than the plotting 
symbol. 
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Fig. 10. — Spectra of continuum components A, B, and E and corresponding spectral indices a, 
where Sv oc i^". The arrow indicates an upper limit on the flux density of component E at 15 GHz. 
The effect of resolution on the spectral index is probably small because the beam dimensions at 
the three frequencies are comparable (Table 1). Flux densities of the continuum components at 
22 GHz are not shown because they were obtained at a different epoch, and the emission is known 
to be time variable. 



reveals variation in the flux density of continuum component A at 22 GHz by a factor of ~ 5 over 
6.2 years. If we compare the results of Sawada-Satoh et al. (2000) and of this study, the variation 
in the flux density of component A at 22 GHz could have been as large as a factor of ~ 52 over 
4.4 years (we note that, although the 22 GHz beam of Sawada-Satoh et al. (2000) is larger, they fail 
to detect component A). In both cases, the timescale of variabiHty is greater than the light travel 
time within the components (< 2 yr). 

Our flux density estimates for component E indicate a spectral index of a = —0.2 between 5 and 
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8 GHz and a spectral index of a < —2.1 between 8 and 15 GHz (Fig. 10). The high-frequency dechne 

in the flux density of component E is too steep to be consistent with optically-thin synchrotron 
emission [a = (1 — p)/2], since it would require p > 5.2, much greater than the normally assumed 
value of p ~ 2.4, where p is the power law index on the electron energy distribution [N{E) oc E~p\. 
However, for a situation in which there is no continuous injection of electrons into the synchrotron 
source and the pitch angle distribution remains conserved, the cooling of high-energy electrons leads 
to a spectral index of a = —{2p-\- 1)/3 above a break frequency ui, = [-6/(1 G)]~^[t/(1 yr)]~^GHz 
(Kardashev 1962). Hence, nominally, a field of ^ 40 mG would give ui, ~ 10 GHz after a cooling 
time of ~ 40 years, the travel time of component E at 0.12 c from the central engine (refer to 
Section 4.1.1). A value oi p > 2.6 would be required to give a spectral index of a < —2.1 above the 
break frequency, consistent with our observations. Using the same magnetic field and extending the 
cooling time by ~ 5.8 years (i.e., the time difference between the two studies), we estimate a break 
frequency at the observing epoch of Middelberg ct al. (2003) to be ~ 8 GHz, consistent with their 
detection of E at 5 GHz and non-dctcction of E at 15 GHz. We note that their spectrum appears 
to be shifted to lower frequencies because it peaks between 1.7 and 5 GHz rather than between 
5 and 8 GHz as in our study. This shift is consistent with an adiabatic expansion of the rapidly 
cooling synchrotron source, for a source lifetime of ~ 40 years, p ~ 2.6, and the ratio of the turnover 
frequencies of ~ 1.4 (corresponding to a shift of roughly 1.5 GHz; note that turnover frequency is 
different from u^; Eq. 13.27 of Kellermann & Verschuur 1988). 

4. Discussion 

In what follows, we propose that the inner parsec in NGC 3079 contains a nearly edge-on, 
massive, thick, and fiared disk. The disk orbits a mass of ~ 2 x 10^ Mq enclosed within 0.4 pc, 
most likely a supermassive black hole, and it is aligned with the kpc-scale molecular disk. The 
disk is most likely self-gravitating, clumpy, and forming stars. Accretion onto the black hole drives 
a jet that is misaligned with the disk rotation axis and is interacting with the dense ambient 
medium. The presence of off-axis aging synchrotron components may suggest that the jet changes 
orientation. The jet may coexist with a pc-scale wide-angle outflow, which is inferred from the 
observation of masers at high latitudes above the disk and which might be related to the known 
kpc-scale superbubble. 

4.1. Analysis of Maser Kinematics 

4.1.1. The Disk Major Axis and the Dynamical Center 

The largely north-south linear distribution of the maser emission, aligned with the kpc-scale 
molecular disk, as well as the segregation of the blue- and the red-shifted emission on the sky are 
suggestive of a nearly edge-on pc-scale molecular disk. A fit to the angular distribution of maser 
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Fig. 11. — Distribution of maser emission (colored circles) and of the 22 GHz continuum (gray 
contours) in the nuclear region of NGC 3079. The color of the maser spots indicates line-of-sight 
velocity in accordance with the bar on the right. The continuum contour levels are —3a, 3a, 5a x 
2^^/-^, where a = 0.14 mJy beam~^ and k = 0,1,2,3,.... The nearly vertical line is the proposed 
disk major axis, while the gray nearly vertical region illustrates the uncertainty in the location and 
orientation of the axis. Also shown in the figure is the outflow trajectory of component A at P.A. 
of 126° lb 5° (as in Fig. 9), where the cone shows the la uncertainty in the orientation of this 
trajectory. The intersection of the outflow direction with the disk major axis yields the location 
for the dynamical center (filled symbol), where 70% and 45% confidence regions for its location 
are illustrated by the central contours. Note that the two easternmost and two westernmost maser 
features (outside the gray area depicting the maser disk) are most likely associated with the outflow. 
The ellipse in the bottom right corner of the figure illustrates beam dimensions. 
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Fig. 12. — Average velocity (v), enclosed total mass M{r), velocity dispersion o", and scale height 
h as a. function of distance from the location of the dynamical center along the disk major axis, 
i.e., impact parameter. At each impact parameter, the four quantities under consideration were 
computed using an averaging aperture of radius 4 mas (0.34 pc) and weighting i*"^ for various values 
of 7, where F is the maser flux density and 7 is denoted by color. The average velocity, velocity 
dispersion, enclosed mass, and disk scale height were computed using {v) = F^] Yli (^i ~ 
Vsys), M{r) = {vfr/G, = Ei-^'^]"^ Ei K - Vsys - {v) f, and h ~ ra/{v), respectively, 
where the sums are over all maser features located within the averaging aperture and r is the 
impact parameter. The curves were computed in a quasi-continuous fashion to emphasize the 
uncertainty in the computed parameters due to discreteness in velocity and spatial sampling. The 
dotted lines show the analytical results for a Keplerian edge-on disk of negligible mass around 
a [0.5,2,4] X 10^ Mq point mass computed using an averaging aperture of 4 mas and under the 
assumption that the maser emission traces completely the midline of the putative disk. 
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emission yields the following disk major axis: a(mas) = (—0.17 it 0.06) x (5(mas) + (0.5 it 0.7) mas, 

where a and b are right ascension and declination, respectively, measured with respect to the 
reference maser feature (Fig. 11). The position angle of the proposed major axis is ~ —10°, in 
close agreement with the position angle of the kpc-scale molecular disk (P.A.~ —11°; Koda et al. 
2002). 

We use the estimated disk major axis and the systematic motion of component A to establish 
the location for the dynamical center. The uncertainty in the location and orientation of the disk 
major axis restricts the position of the dynamical center to lie within the gray nearly vertical band 
in Fig.ll. The systematic motion of component A (Fig.9), on the other hand, limits the location of 
the dynamical center to P.A.^^ —54° it 5° as measured with respect to A and illustrated as a cone in 
Fig.ll. The intersection of these two loci yields (a, 5) = (1.9 it 1, —8.2 it 2) for the mean position of 
the dynamical center (shown in Fig. 11). Note that the proposed position of the dynamical center 
lies in close proximity to the intersection of the disk major axis with the line joining components A 
and B, which was expected since the trajectory of A is nearly parallel to the line joining A and B. 
The position of the dynamical center proposed in this study is similar to that proposed in Trotter 
et al. (1998) but significantly different from that proposed by Sawada-Satoh et al. (2000) (Fig. 1). 



J^.1.2. Disk Structure 

To study the disk structure traced by the maser emission, we computed the average velocity 
(f), the enclosed total mass M(r), the velocity dispersion cr, and the disk scale height ft. as a 
function of distance r from the putative dynamical center along the major axis of the disk, i.e., 
the impact parameter (Fig. 12). In this computation and the following analysis, we explicitly 
assume that the maser emission traces the kinematics of the underlying molecular material, as is 
indeed the case in the archetypal maser galaxy NGC 4258. To compute the four average quantities 
at each impact parameter, we used an averaging aperture of radius 4 mas (0.34 pc) and weighting 
function ^ where F is the maser flux density and < 7 < 1. This particular weighting scheme 
was adopted to demonstrate that the general conclusions inferred from this analysis are relatively 
insensitive to the observed broad range of the maser flux density (Fig. 3). Moreover, at each 
impact parameter, the average velocity and velocity dispersion were computed in accordance with 

{v) = [Ei FiV^ Ej Fi (^i - Vsys) and cr^ = F^]'^ Y.-, ^7 - ^'«y*^ - (^'))^ respectively, where 
the sums are over all maser features located within the averaging aperture. Assuming Keplerian 
rotation and spherical symmetry, we estimated the enclosed mass as M(r) = {v)'^r/G, where r is the 
impact parameter. Under the assumption of hydrostatic equilibrium and axisymmetric potential, 
the disk scale height is given by /i ~ ra/{v). Note that the two easternmost and two westernmost 
maser features have been omitted from this analysis because they do not appear to be associated 
with the disk (Fig. 11). Their inclusion would not change our conclusions significantly. 
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4.1.3. Rotation, the Enclosed Mass, and the Disk Mass 

The computed average velocity displays a red-blue asymmetry about the adopted Vsys and 
estimated dynamical center and is thus consistent with rotation (Fig. 12). The Keplerian rotation 
law yields a mass of ~ 2 x 10^ enclosed within 0.4 pc (5 mas). Trotter et al. (1998) reported 
a similar enclosed mass, ~ 10^ Mq, from consideration of the blue-shifted emission only. In our 
case, however, both the blue- and the red-shifted emission provide independent and consistent 
estimates of enclosed mass. We note that since the maser features may lie somewhat away from 
the disk midline or the disk might be somewhat tilted from edge-on, the computed enclosed mass 
is a lower limit, but is probably correct to within factors of order unity. For a central mass of 
2 X 10^ Mq and an inner radius of 0.4 pc, the maximum centripetal acceleration (i.e., velocity drift) 
and the maximum proper motion would be 0.05 kms~^ yr"^ and 2 fiasyv~^, respectively. The much 
larger apparent velocity drifts and motions (Fig. 5) are most likely due to variability of separate 
gas clumps and macroscopic random motions among clumps rather than global kinematics, which 
makes a measurement of velocity drifts and proper motions due to the gravitational potential of a 
central mass difficult if not impossible. 

Since the rotation curve of both blue- and red-shifted features is flat, i.e., M(r) oc r, the 
mass of the circumnuclear pc-scale disk is significant with respect to the central mass. In fact, the 
disk mass within 0.7 pc, the distance from the dynamical center to the brightest maser feature, is 
roughly Mrf(0.7 pc) ~ M(0.7 pc) - M(0.4 pc) ~ 10^ Mq. The mass of the entire disk traced by the 
maser emission (< 1.3 pc) might be as large as 7 x 10® Mq, but this is uncertain by a factor of a 
few. Interestingly, if we extrapolate M(r) oc r to larger radii, we estimate a mass of ~ 4 x 10^ Mq 
enclosed within 76 pc, which is in surprisingly good agreement with the 7 x 10^ Mq dynamical mass 
inferred from CO observations (Koda et al. 2002). 

4.1.4. Mean Mass Density, Eddington Ratio, and Mass Accretion Rate 

The mean mass density corresponding to 2 x 10^ enclosed within 0.4 pc is 10^-^ Mq pc""^. 
The relatively high mean mass density for NGC 3079 is suggestive of a massive central black hole as 
opposed to a dense star cluster (e.g., Maoz 1995a), which is consistent with the X-ray observations 
of the nucleus. The estimated enclosed mass of ~ 2 x 10® Mq is in agreement with proposed correla- 
tions between bulge velocity dispersion and black hole mass (Gebhardt et al. 2000a, 2000b; Ferrarese 
& Merritt 2000; Ferrarese et al. 2001). If we adopt Mbh = 1-2 x 10^ Mq [a/ {200 km s"^)]^-^^ (Geb- 
hardt et al. 2000a, 2000b), then the velocity dispersion of the bulge, 40kms""^ < a < 160kms~^ 
(Shaw, Wilkinson, & Carter 1993), is consistent with a black hole mass of 10^-5-^-^Mo, although 
the heavy dust obscuration might bias the measured dispersion. 

To estimate the mass accretion rate as well as the accretion timescale, we first approximate 
the luminosity of the nucleus. Since the 2 — lOkeV luminosity is probably 5% of the AGN 
bolometric luminosity (e.g., Kuraszkiewicz et al. 2003, Elvis et al. 1994), we obtain an AGN 
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bolometric luminosity of ~ 5 x 10^~^° Lq from the X-ray data analysis reported by lyomoto et al. 
(2001). The Eddington luminosity of a 2 x 10^ Mq object is 7 x 10^° Lq and, assuming that all of 
the enclosed mass is concentrated in a supermassive black hole, the approximate luminosity of the 
central engine yields an Eddington ratio of 0.08 — 0.8, which is consistent with the 0.01 — 1 range 
obtained for Seyfert 1 galaxies, representative supermassive black hole systems (e.g., Padovani 1989; 
Wandel 1999). Assuming a standard accretion efficiency of ~ 0.1 (Prank, King, & Raine 2002; see 
also Marconi et al. 2004), we obtain a mass accretion rate of M = 0.007 LboI,io -^©year"^, where 
LboI.io is the AGN bolometric luminosity in units of 10^'^ Lq. The disk mass computed above yields 
an average accretion timescale oi t = M^{r)/M ~ 10* ^BoI 10^®^'-^ ^ radius of 0.7pc (Table 4). 



Table 4: Parameters of the central engine and pc-scale disk in NGC 3079. 



Quantity 



Estimate 



Central mass (r < 0.4 pc) 

Mean mass density 

Disk mass (0.4 pc < r < 0.7 pc) 

Disk mass (0.4 pc < r < 1.3 pc) 

Disk scale height 

AGN bolometric luminosity 

Eddington ratio 

Accretion rate'' 

Accretion timescale''''' 

Toomre Q-parameter 

Clump size* 

Clump mass" 

Jeans mass 

Rotation period'^ 

Clump cooling timescale" 

Clump free-fall timescale 

Clump collision timescale 

Roche limit" 



O.i 



~ 2 X 10" Mq 
~ 10^-^ Mq pc"^ 
~ 10® Mq 

< 7 X 10® Mq 

0.05 pc < ft < 0.5 pc 
- 5 X 108-1° Lq 
0.08-0.8 
007 LboI. 10 Mq year-i 
lO^-^-BoLio years 
0.01 <Q < 0.02 

< 0.006 pc 

< 5 X 10^ Mq 
0.3 - 53 Mo 

4 X 10^ years 

< 60 years 

W*-^ years 
n > 5.3 X 10* cm"^ 



"Computed at the location of the brightest maser feature, r = 0.7 pc. 
''LboI,w is the AGN bolometric luminosity in units of 10^° Lq. 



4.1.5. Disordered, Thick, and Flared Disk 

The disk in NGC 3079 is different from the archetypal Keplerian disk in NGC 4258 in that 
the velocity structure of the former is much more disordered. The rotation traced by the maser 
emission in NGC 3079 is characterized by relatively large velocity differences across relatively 
small areas on the sky. For instance, the velocity dispersion (a) as computed using an averaging 
radius of 4 mas (0.34 pc) ranges from 20 to 80kms~^ (Fig. 12). The maximum velocity difference 
between neighboring maser features is 144kms~^ across a region as small as 0.1 pc. The large 
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velocity dispersion is most likely indicative of macroscopic random motions among the molecular 
clumps responsible for the maser emission rather than turbulence within the clumps (Section 4.2). 
As already noted above, the variability in the angular distribution and in the spectrum of the 
maser (Figs. 3, 5) is a direct consequence of these macroscopic random motions. Although a 
significant fraction of the orbital velocity, the computed velocity dispersion is everywhere much 
smaller than the escape velocity of a ~ 2 x 10^ central object {vescape > llOkms"^ in the 
region supporting maser emission) and thus has little impact on the stability of the underlying 
rotating structure. Although we favor the macroscopic random motions of clumps as the origin of 
the observed high dispersion, we note the small size of the dominant maser feature (< 5 x 10"*^^ cm) 
for which a ~ 14kms^^, and we speculate that there may be some regions in which turbulence is 
significant but does not interfere with maser action (e.g., Wallin, Watson, Sz Wyld 1998). 

The computed velocity dispersion can be used to infer the disk thickness. Because cr is a 
considerable fraction of the orbital velocity, the rotating structure is most likely geometrically thick 
(e.g., Alves Sz Nelson 2000; Wainscoat, Freeman, &; Hyland 1989). Moreover, since a does not 
decrease with the impact parameter, the disk scale height ought to increase with the distance 
from the center; i.e., the disk is probably flared. In particular, the scale height, computed using 
h ~ ra/{v) appropriate in the case of hydrostatic equilibrium, is on the order of 0.05 pc at an 
impact parameter of 0.5 pc and might be as large as /i ~ 0.5 pc at a radius of 0.9 pc (Fig. 12). The 
computed scale height is in reasonable agreement with the apparent thickness of the disk inferred 
from the dispersion of the maser spots about the disk major axis, h ~ 0.25 pc at r ~ 0.7 pc (see 
Fig. 11). 



4.2. Clumpy Star-Forming Disk 

The Toomre-Q parameter characterizes the stability of the accretion disk, and we can estimate 
it based on our measurements and an assumption about the temperature of the maser medium. 
The Q parameter is given by Q = VLCg/ {ttGTi)^ where VL is the angular rotation rate, S is the 
surface density, and Cg is the sound speed. If we assume that the mass distribution in the accretion 
disk has circular symmetry, then S = {dM / dr) / {2TTr) and we can write Q = 2vCs / {G dM / dr) . 
Using M(r) = v^r/G and the fact that the orbital velocity of u ~ 150kms~^ is approximately 
constant with radius (Fig. 12), we estimate a radius independent value of roughly dM/dr = 
/G 3 X 10^"*^ gcm~^. The sound speed in a neutral medium is given by cs = 0.04Ti/2ijms-^ 
(e.g., Maoz & McKee 1998), and most models of maser emission from water vapor require the gas 
temperature to lie in the range 300 — 1000 K (e.g., Desch et al. 1998). With these parameters we 
obtain Q = 0.01 — 0.02, where the range reflects only the uncertainty in temperature. Except for 
possible temperature variation, our estimate of Q is independent of radius since both the orbital 
velocity and dM/dr are to first order constant with radius. Since Q is significantly less than unity, 
the disk appears to be gravitationally unstable in the region supporting the maser emission. The 
flat rotation curve, which implies signiflcant accretion disk mass with respect to the black hole 
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mass, is fundamentally responsible for the low value of Q. The value of Q we estimate is rather 
robust if our fundamental assumption that the velocities of the masers are due to the gravity of 
the enclosed mass holds true. 

From our estimate of dM/dr, we parameterize the disk surface density as E ~ 250 (r/0.7 pc)~^ gem' 
i.e., about 430 and 130gcm~^ at the inner and outer radii of 0.4 and 1.3 pc, respectively. (The 
empirically derived power law index of —1 is not unreasonable; for instance, E oc r~^/^ for a ra- 
diatively cooled, gas pressure-dominated cc-disk with a constant rate of accretion throughout the 
disk.) Using the scale height in Fig. 12, our estimate for disk surface density, and S = 2h{p), 
we compute mean density that varies from about 0.3 x 10^ to 3 x 10^ mn cm~^ over the observed 
disk. The densities deduced here lie comfortably within the normally required range of H2 number 
densities of lO'''"^'' cm^^ for maser emission (e.g., Desch et al. 1998). Our determination here of 
the average density of a maser medium is probably the first accurate determination of this quality 
for water in general. 

The ultimate fate of an accretion disk characterized by Q < 1 depends on the balance between 
cooling and heating. It has been argued that the disk will maintain Q ~ 1 by rearrangement of sur- 
face density to reduce S or by heating through turbulence (Lin &; Pringle 1987; Hure 2000; Gammie 
2001). However, the energy required to maintain Q ~ 1 is prohibitive (Goodman 2003) and, in the 
case of NGC 3079, would require temperatures as high as 10^ K. Furthermore, cooling timescales 
might be sufficiently short that a Q 1 criterion is unsupportable (Shlosman &: Begclman 1989; 
Monaghan & Lattanzio 1991), in which case the disk would fragment into clumps (e.g., Gammie 
2001; Kumar 1999). 

We consider the clumpy disk model in detail by computing clump size and mass at a repre- 
sentative location within the disk, the distance from the dynamical center to the brightest maser 
feature, 0.7 pc. The gas number density required for gravitational collapse in a strong tidal field 
IS given by the Roche limit, n > 5.3 x 10^ (M(r)/2 x 10^ M©) (0.7 pc/rf cm , which is consis- 
tent with the H2 number density required by the presence of the maser emission. Furthermore, as 
derived by VoUmcr, Beckcrt, & Duschl (2004), an upper limit on the radius of a self-gravitating 
clump subject to tidal shear is Rc < TrCg/V^O., from which we obtain Rc < 0.006 pc and an upper 
limit on the clump mass of 5 x 10^ M©. It is instructive to compare this upper limit to the Jeans 
mass, mj = 

(7rc2/G)3/2p-V2^ which is 0.3- 53 M© given the conditions necessary to support maser 
emission. Since the Jeans mass and the Roche limit are consistent with the upper limit on cloud 
mass and with H2 number density inferred from the presence of the maser emission, respectively, 
fragmentation and star formation are expected to occur within the clumps. 

In addition to the Jeans mass analysis, straightforward considerations of energetics also sug- 
gest that individual clumps are susceptible to collapse. The time in which a clump will radi- 
ate its binding energy via black-body emission is given by tc ~ (GM^ / Rc)/ {4:nR^aT^) 
(GMcmH(nH2))/(3(Tr^) where Mc, Rc, Tc, and (nna) are the clump's mass, radius, effective tem- 
perature, and average H2 number density, respectively (e.g., Goodman 2003). In our case, Mc < 
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4 X 10^ M0, Tc > 300 K, and {nn^} < 10^° cm-^, which yields tc < 60 years. However, if Une emis- 
sion is the dominant coohng mechanism, then the timescale becomes tc ~ {G / Rc) / {4tt R^A/S) ~ 
(2.3)^ X (47rGm|^i?^(nH2)/3)(A/(nH2))~^, where A/{nu2) total cooling rate per H2 molecule 

in ergss"^. Using the cooling rates in Neufeld et al. (1995) and Rc < 0.006 pc, the expression above 
is maximized for {71112) — 10^*^ cm~^ and = 300 K, which results in tc < 40 years. If cooling of 
the gas by cold dust grains dominates, then, using Rc < 0.006 pc, dust temperature of ~ 200 K 
(Collison Sz Watson 1995), and the cooling rate per Tuiit volume of HoUenbach & McKee (1989), 
we obtain a cooling timescale of at most 5 years. Thus, the time in which a clump will radiate all 
of its binding energy via either black-body radiation, line emission, or heating of cold dust is much 
less than the accretion timescale and even the disk rotation period (Table 4); consequently, clumps 
collapse on their free-fall timescale, resulting in star formation. In our case, the free-fall timescale 
is ^/Rf/GM Ki (37rGmijnH2)~°'^ = 10^-^~^ ° years, depends only on the hydrogen number density 
implied by the conditions necessary for maser amplification, and is four orders of magnitude smaller 
than the computed accretion timescale. 

Star formation can also be triggered by collisions among clumps. W49N, the most luminous 
water maser in our Galaxy, provides probably the best example of how clump-clump interaction can 
induce 0-type star formation (Serabyn, Giisten, &: Schulz 1993). To estimate the collision timescale, 
we must obtain the mean clump separation. Using the scale height in Fig. 12, we estimate the 
volume of the entire disk as traced by the maser emission to be roughly 2pc^. Mass of the entire 
maser disk of ~ 7 x 10^ Mq and a 5 x 10^ Mq upper limit on clump mass yield > 1.5 x 10^ clumps in 
the disk and mean clump separation of Is < 0.03 pc. Hence, the collision timescale, ~ {Ig/Rc)^ 
(Kumar 1999), is ^ 10"^^^ years and in fact depends only on disk size, disk mass, and H2 number 
density inferred from the presence of the maser emission. Since the collision timescale is much 
less than the accretion timescale, collisions among clumps might also be a significant process in 
triggering star formation and overall evolution of the disk. 

In summary, since Q <C 1, the pc-scale circumnuclear disk in NGC 3079 is most likely gravita- 
tionally unstable and therefore clumpy. Evidence for a self-gravitating disk of mass comparable to 
the central mass has also been found in NGC 1068 (M^ Mbh = 8.0 x 10^ M©; Lodato & Bertin 
2003), although NGC 3079 provides a more extreme example of the phenomenon (M^ 7x 10^ Mq 
vs. Mbh ~ 2 x 10® Mq). Furthermore, the Jeans mass, the Roche limit, as well as cooling and 
collision timescales are suggestive of star formation within the clumps through either clump collapse 
or clump collisions. Once a clump collapses, star formation proceeds via accretion of gas from the 
pc-scale accretion disk onto the Hill sphere of the stellar object, as considered in detail by Milosavl- 
jevic & Loeb (2004). The star formation efficiency of giant molecular clouds (M = 10^~^ Mq) and 
dark molecular cores (M = 1 Mq) ranges from 0.1 % to 5 % while the star formation timescale and 
the lifetime of hot massive stars are on the order of 10® years (e.g.. Blitz & Thaddeus 1980; Larson 
1982; Wilking & Lada 1983; Myers 1985; Silk 1985, 1987). Thus, the long accretion timescale 
and the large potential clump mass computed above are consistent with the relatively slow and 
inefficient process of star formation (Table 4). If the star formation were to exhaust the available 
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gas (which, in addition to the star formation rate, also depends on the presence of any gas feeding 

mechanisms), then the massive accretion disk in NGC 3079 would transform into a stellar disk, a 
transition that has been proposed to explain the presence of young stars in the immediate vicinity 
of the black hole in the Galactic center (Levin & Beloborodov 2003; Milosavljevic &: Loeb 2004; 
Nayakshin & Cuadra 2004). 

4.3. Pumping of Maser Emission by the Central Engine 

Since the nucleus of NGC 3079 contains a compact hard X-ray source, irradiation of molecular 
gas is a plausible means of exciting maser emission (e.g., Neufeld, Maloney, Sz Conger 1994). In 
maser sources such as NGC 4258 and Circinus, irradiation over a wide range of radii is achievable 
because the disks are warped (Herrnstein, Greenhill, & Moran 1996; Greenhill et al. 2003). In 
NGC 3079, the disk in which the masers lie does not appear to be warped, in which case maser 
excitation over the observed range of radii (0.4 — 1.3 pc) depends on penetration of X-rays in the 
disk plane. However, because maser action requires high densities, the disk must be inhomogeneous, 
which is consistent with our earlier stability arguments. For a minimum H2 density of lO^cm"^, a 
column within the disk becomes Compton thick for lengths greater than ~ 0.02 pc. However, for a 
line-of-sight filling factor of ^ 2%, the column could pass X-rays to a radius of ^ 1.3 pc. 

For an irradiated slab of gas, Collison &: Watson (1995) estimate a maser emission rate of 
~ 4000 Lq per pc^ of surface area beamed in the plane of the slab. For a model spherical gas clump 
of radius ~ 0.005 pc that is irradiated on one side, the integrated output is on the order of 0.3 Lq, 
which is much less than the isotropic luminosity inferred from the strength of individual maser 
features (e.g., 131 L©; Fig.3). We suggest instead that the emission from the observed masers is 
narrowly beamed along our line of sight by the overlap of clumps with similar Doppler velocities, 
whereby the integrated output of each is directed toward tis (e.g., Deguchi &; Watson 1989; Kartje 
et al. 1999). To support the peak observed isotropic luminosity, we require two clumps, each 
0.005 pc in size, separated by on the order of 0.1 pc, which is reasonable in light of the ~ 1.3 pc 
radius of the disk. We note that amplification of background continuum emission by individual 
clumps can also generate strongly forward beamed maser emission, but we do not detect continuum 
emission in the vicinity of the disk above 0.84 mJy (6cr) at 22 GHz. On the other hand, the two 
easternmost and westernmost maser features do lie in close vicinity to continuum components B 
and C, and we speculate that maser emission observed away from the disk, where the density of 
clumps is probably greatly reduced, may be the result of continuum amplification. 

Under the assumption that clumps are optically thick, irradiation of individual clumps by a 
central source creates a dissociation region on the inward facing side, which results in anisotropic 
emission of microwave photons preferentially along the length of the dissociation region. This 
pattern may be used to explain why we observe maser emission in the vicinity of the disk midline 
but not close to the line of sight to the central engine. Clumps on the front and back sides of the disk 
beam maser photons along axes that are substantially inclined to our line of sight. Furthermore, 
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we speculate that the observed variabiUty of maser spectrum (Nakai et al. 1995; Baan & Haschick 
1996; Hagiwara et al. 2002) and the observed variation in angular structure (Fig. 5) may be the 
result of pairs of clumps moving in and out of alignment along the line of sight. The presence 
of the maser feature at ~ 956 km in the spectrum of NGC 3079 since the discovery of the 
maser (Henkel et al. 1984; Haschick Sz Baan 1985; Baan Sz Haschick 1996) and the a ~ 30kms~^ 
dispersion due to bulk random motion (Fig. 12) place a lower limit on the clump size of ~ 0.001 pc, 
consistent with the upper limit (Table 4). Year-to-year variation in the flux density of some maser 
features might be due to incomplete overlaps on the sky of the X-ray dissociation regions. 

4.4. Outflow 

4.4- 1- Components E and F 

Component E, identified by Kondratko et al. (2000), and component F, identified by Middel- 
berg et al. (2003), might be remnants of jet emission along former axes of a wobbling jet, currently 
traced by components A, B, and C. The apparent decrease in the break frequency of components E 
and F with distance from the central engine (Middelberg et al. 2003) and the shift of the spectrum 
of component E toward lower frequencies with time (Section 3.2) are indirectly supportive of the 
hypothesis that both components are rapidly cooling and adiabatically expanding remnants of jet 
emission. Although the structure of the continuum is difficult to understand in detail, we speculate 
that components A and B are regions where the jet, in its current orientation, interacts with a 
dense ambient medium, an inference supported by their GPS-like spectra (Section 3.2). 

Alternatively, considering the location of E (Fig. 13) and of F (Middelberg et al. 2003) within 
the opening angle of the kpc-scale superbubblc, these two components might provide evidence for a 
pc-scale wide-angle outflow, most likely an inward extension of the kpc-scale superbubble. In this 
interpretation, the two VLBI components could be moving knots or shocks as the ionized wide-angle 
wind interacts with the dense ambient medium along the edges of the flow as seen in projection. 
If this hypothesis is true, then new continuum components are expected to appear on both sides 
of the central engine but within the opening angle of the superbubble. However, it is currently 
not possible to distinguish between the two models proposed above, and more sensitive images 
and proper motion studies of the continuum components extending over a longer time frame are 
necessary. Given the possible source lifetimes and transverse speeds (Section 3.2), time baselines 
on the order of a decade are necessary. 

On the other hand, it is unlikely that E and F arc radio supernovae associated with the nuclear 
starburst. The spectral index of E, a < —2.1 between 8 and 15 GHz (Fig. 10), is much steeper 
than what has been observed in radio supernovae (a > —0.9; Weiler et al. 1986; Weiler, Panagia, 
& Sramek 1990; Allen & Kronberg 1998; McDonald et al. 2001, 2002; Bartel et al. 2002) . The 
extremely steep spectrum of F, a ~ —6.1 between 1.7 and 2.3 GHz, is also not compatible with a 
radio supernova. Furthermore, the radio supernovae in the prominent starburst galaxy M82 scaled 
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to the distance of NGC 3079 would have flux densities (< 0.1 mJy) below our detection limits at 5 
and 15 GHz (McDonald et al. 2002). 

4-4-^- Wide-Angle Outflow on Parsec Scales 

The existence of a wide-angle outflow in NGC 3079 was hypothesized in previous studies of 
the nucleus. Based on VLA data, Duric & Seaquist (1988) first suggested that a wide-angle wind 
is responsible for a bipolar "figure eight" structure along the minor axis of the galaxy in radio 
continuum. A wide-angle outflow on pc-scales has been suggested to explain a blue-shifted OH 
absorption component and two weak OH emission features detected from the nucleus (Baan & 
Irwin 1995; Hagiwara et al. 2004). By considering the energetics and morphology of the kpc-scale 
superbubble observed in Hq with the Hubble Space Telescope (HST), Cecil et al. (2001) argued 
that the bubble is inflated by a wide-angle outflow rather than a processing jet. Wide-angle outflows 
seem to be quite common among Seyfert galaxies (Colbert et al. 1996a,b, 1998), and one has been 
imaged on pc-scales in the Circinus galaxy, the nearest Seyfert 2 nucleus (Greenhill et al. 2003). 
Guided by this result, by the presence of molecular gas as dense as the disk gas at high latitudes 
above the disk (i.e., the two easternmost and two westernmost maser features in Fig. 11), and by 
the existence of a kpc-scale superwind, we suggest a wide-angle outflow on pc-scales in NGC 3079. 

It has been suggested that winds are driven by photoionized evaporation of matter from an 
inner surface of a torus (Krolik & Kriss 2001), by radiation and gas pressure acting on an accretion 
disk (Murray et al. 1995), or by the magneto-centrifugal uplift of gas and dust from an accretion 
disk (Emmering, Blandford, &: Shlosman 1992; Konigl &; Kartje 1994, Kartje et al. 1999). In 
addition, ram pressure of a wide-angle outflow can entrain clumps at the surface of an accretion 
disk. If the wind is sufficiently dense to effectively shield the clumps from the central engine and the 
clumps can be confined, then they can potentially rise to high latitudes above the rotating structure 
(see Fig. 7 in Kartje et al. 1999) while still maintaining, to some degree, the rotational velocity 
imprinted by the parent disk. In light of this, the fact that the line-of-sight velocities of the two 
easternmost and the two westernmost features (Fig. 11) reflect the velocity of the most proximate 
side of the disk can be explained if these clumps were uplifted from the disk surface and carried 
to high latitudes by a dense wide-angle outflow. We thus consider the four maser features located 
significantly out of the plane of the disk as indirect evidence for a pc-scale wide-angle outflow, likely 
an inward extension of the kpc-scale superbubble. The fact that some clouds attain large latitudes 
above the disk surface while others remain within or in a close proximity to the thick disk might be 
due to a mixture of heavy and light clumps within the disk. We note that, due to beaming effects 
of maser emission, the actual number of the clumps uplifted might be much greater than the four 
clumps that are observed to be associated with the wind. 

Although the origin of the putative pc-scale wide-angle outflow is unclear, we suggest that it 
is unlikely to be driven by star formation in the pc-scale accretion disk. The mechanical luminos- 
ity required to blow out a superbubble from the galactic disk in NGC 3079 is ~ 3 x 10^^ ergs s~^ 
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(Veilleux et al. 1994). Following Strickland et al. (2004a), such a luminosity would require ~ 
10^ (t/Myr) supernovae corresponding to a cluster of mass 5 x 10^ (t/Myr)MQ, assuming a ther- 
malization efficiency of ~ 1, an energy release of 10^-*^ ergs per supernova, and a starburst lifetime t 
equal to the dynamical age of the superbubble, ~ 10^ years (Veilleux et al. 1994). The predicted 
cluster mass is an order of magnitude greater than the 21, 000 — 79, 000 Mq mass of one of the most 
massive star clusters known, R136 in 30 Doradus, with a half-mass radius of 0.5 pc (Campbell et 
al. 1992; Brandl et al. 1996; Massey &: Hunter 1998). Such massive clusters are thus extremely 
rare, constituting only a fraction of ~ 5 x 10^^ of all star clusters (Strickland et al. 2004a). More 
reasonable values of the starburst lifetime (~ 40 Myr) would require an even more extreme cluster 
mass and number of supernovae, effectively excluding a pc-scale cluster as a progenitor of the kpc- 
scale superbubble. Supernovae localized in multiple starburst regions and distributed over larger 
scales (e.g., kpc) might still contribute significantly to the wide-angle outflow, although the degree 
of this contribution is uncertain (Strickland et al. 2004b,a). 

One consequence of the broad outflow that we hypothesize — or the canted jet marked by 
continuum components A, B, and C — may be the uniform weakness of maser emission south 
of the dynamical center. Since AGN outflows are ionized, differential attenuation of the maser 
emission is readily achieved for free-free absorption with a small difference in electron density and 
corresponding emission measure. If the disk is tilted slightly from edge-on such that the lines of sight 
to the maser emission traverse outside of the disk (Fig. 13), the observed two order of magnitude 
flux density difference between blue- and red-shifted features may be explained by a factor of at 
most 2 difference in local electron number density (provided that r > 1 in the direction of 
the blue-shifted maser emission). In the case of a canted jet, an enhancement in density may be 
readily achieved given that the jet is inclined toward the apparently weaker masers (see also Trotter 
et al. 1998), but a density inhomogeneity in the wide-angle outflow might also be the cause of the 
differential attenuation. However, in either case, it is interesting to note that the balance of blue- 
and red-shifted emission can be readily governed by factors extrinsic to the maser emitting regions 
in the case of both NGC3079 and NGC4258 (Herrnstein et al. 1996). This may be generalizable 
to other maser sources and contrasts with the models of Maoz (1995b) and Maoz & McKee (1998). 



5. Summary 

Water maser emission in the active nucleus of NGC 3079 has recently been recognized to cover 
a ~ 450kms~^ range centered on the systemic velocity of the galaxy. We have mapped for the first 
time maser emission over this entire velocity interval. We have also imaged non-thermal continuum 
emission between 5 and 22 GHz that arises in close proximity to the molecular gas underlying the 
masers. Based on the analysis of the spectral-line and continuum maps, we conclude the following: 

1. The largely north-south distribution of maser emission, aligned with a known kpc-scale molec- 
ular disk, and the segregation of blue- and red-shifted maser emission on the sky strongly 
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Fig. 13. — Left: [Nil] + Ha image of the superbubble (filled color contours) obtained with HST 
(Cecil et al. 2001) and the soft X-ray image of the same (black contours) obtained with Chan- 
dra. The X-ray image has been convolved with a Gaussian of 0.64" FWHM. The cross indi- 
cates both the VLBI-derived position of the maser source with respect to the X-ray emission 
and the 0.6" uncertainty in the Chandra absolute astrometry (Chandra X-Ray Center Web site: 
http://cxc.harvard.edu/cal/ASPECT/celmon/). Note that the position of a nearby soft X-ray 
emission peak, aaooo = lO'^Ol^ST^Ql ± 0.06, (^2000 = 55°40'47'.'7 ± 0.6, is offset from the dynamical 
center of the maser disk by ~ 1". The alignment of the two images has been accomplished by 
computing their two-dimensional cross-correlation. Although the uncertainty in the HST absolute 
positions might be as large as a few arcseconds, Monte Carlo simulations of the two-dimensional 
correlation, performed by repeatedly adding Gaussian or Poisson noise to the original images, in- 
dicate that the two images are registered to within 0.16". Right: Proposed model for the nuclear 
region of NGC 3079, a VLBI image of 8 GHz continuum (see Fig. 7), and a map of 22 GHz maser 
emission (shown as circles color-coded by Doppler velocity). We propose a thick and flared disk 
(shown in a cross-section) that is slightly tilted with the west side being closer and that is aligned 
with both the kpc-scale molecular disk and the axis of the kpc-scale superbubble. The color-coding 
indicates the Doppler shift of the disk material. The water maser emission traces the disk and may 
arise from X-ray irradiated clumps within the disk, which is unstable to fragmentation. An off-axis 
jet traced by the continuum components A, B, and C probably coexists with a wide-angle outflow 
(yellow), most likely an inward extension of the kpc-scale superbubble. The outflow may result in 
an uplift of clumps from the disk surface, which would explain the high-latitude maser features and 
their association in velocity with the most proximate side of the disk. 
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support the model of Trotter et al. (1998) in which the masers trace a nearly edge-on molecular 
disk about 1 pc in radius. 

2. The dynamical mass enclosed within a 0.4 pc radius is 2 x 10^ Mq, and for a 40 to 160kms~^ 
stellar bulge velocity dispersion is consistent with correlations between central mass and 
dispersion reported for broad samples of galaxies (Gebhardt ct al. 2000a, 2000b; Fcrrarese 
& Merritt 2000; Ferrarese et al. 2001). The ratio of bolometric to Eddington luminosity is 
0.08 — 0.8, consistent with accretion efficiencies of 0.01 — 1 for Seyfert 1 systems. 

3. The disk rotation curve is relatively flat, which is consistent with a mass of < 7 x 10^ Mq 
between radii of 0.4 and 1.3 pc. 

4. The angular distribution of maser emission is not as well ordered in NGC 3079 as in other 
"maser galaxies" (e.g., NGC4258). The velocities of adjacent clumps of maser emission can 
differ by tens of km s~^ on scales of ~ 0.1 pc. As a result, we suggest that the disk is relatively 
thick {h/r ~ 0.1 to 0.5) and may be flared. 

5. Based on stability, cooling, and timescale arguments, we argue that the disk is self-gravitating, 
clumpy, and appears to meet the necessary conditions for star-formation. The maser emission 
most likely occurs as a result of clumps irradiated by X-rays from the central engine and 
overlapping on the sky, yielding long gain paths and narrow beam angles. 

6. We report detection of a very steep spectrum synchrotron component (E) that is not coUinear 
with the previously claimed compact jet. The spectrum is consistent with an aging electron en- 
ergy distribution. The observation of a "relic" favors the hypothesis that the jet in NGC 3079 
has changed direction. In addition to the jet, we observe molecular gas that is dense enough 
to support maser emission at high latitudes above the disk. From this, we infer that the jet 
coexists with a wide-angle outflow originating at parsec or smaller scales. 

The proposed model can be tested through further VLBI study of the NGC 3079 AGN. More 
sensitive images and proper motion studies of the continuum components, with a longer time 
baseline, may help to determine the nature of components E and F. If they are shocks in a wide- 
angle outflow rather than remnants of jet emission, then new continuum components are expected 
to appear within the arc subtended by the superbubble. The presence of star-formation in the 
inner parsec can be confirmed with an ultrasensitive search for M82-like radio supernova remnants. 
A proper motion and monitoring study of the maser emission would be challenging (i.e., motions 
< 2 /xasyr"^) but could be used to corroborate the proposed geometry for the inner parsec, although 
it might be difficult to disentangle motions due to kinematics from apparent motions due to local 
effects, such as clump alignments. Hosting a central engine, a jet, a thick and self-gravitating 
accretion disk, and possibly a wide-angle outflow and star-formation, the nucleus of NGC 3079 
constitutes a nearby laboratory for diverse astrophysical phenomena and is a strong candidate for 
the study of the starburst-AGN connection and of the complex interactions between gas, stars, and 
supermassive black holes. 
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